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Abstract — This article investigates the performance of the generator excitation loop equipped with a PID controlled AVR 
system, and of the dam automaton consisting of a PID controlled water flow regulation system in the southern interconnected 
grid (SIG) of the Cameroonian power system, when subjected to small perturbations, using linearized state-space models of 
the power plant. A fifth order model of the synchronous generator with a high gain excitation system and AVR is considered. 
A detailed formulation of equations comprising the mechanical and electrical swing dynamics of the turbine-generator unit 
and the load is realized; a multivariable state-space non linear model of the one area single machine system model is 
obtained. For small-signal performance analysis, the model is linearized around an operating point. The article also 
presents a state-space model of the water flow system. The effects of PID controlled AVR and water flow regulation systems 
are examined through performing of extensive MATLAB simulations to analyze the behavior of the proposed models 
following small disturbances. The simulation results presented in this paper are obtained using a MATLAB computer 
program developped by the authors; they provide useful insight into the dynamic behavior of the Cameroonian SIG, 
including stability, speed of response and steady-state accuracy. The paper establishes that the classical control systems do 
not have very good performance: long settling time, high overshoot, relativly slow response, and many damped oscillations. 

Keywords — AVR, excitation system, PID controller, state-space modeling, small-signal stability. 

I. Introduction 

Stability is one of the important issues of safe power system operation. Repeated global power outages due to power system 
instability reveal the importance of the issue [1]. Power system oscillations were first observed as soon as two or more 
synchronous generators were connected together to provide more generation capacity and reliability. Originally, the fairly 
closely connected synchronous generators were observed to swing against each other at low frequency; these spontaneous 
oscillations appeared in mechanical variables like rotor angle and speed [2], Since synchronous generators are rotating 
electromechanical devices, mechanical oscillations of their rotor were transmitted by electromagnetic induction to electrical 
variables like bus voltages, line currents and power. The oscillations are initiated by variations in generation and custom 
loads, which act upon the systems as perturbations or disturbances. Because power systems have an almost continuously time 
varying nature - load and thus power demand vary in time, they always suffer from oscillations. Power system oscillations 
are therefore inevitable and are a characteristic of the system. The perturbations could be normal or abnormal. Normal 
perturbations are often of small magnitude, and abnormal perturbations like loss of a large generator or custom load, and 
short-circuit on a transmission line or in a substation are of large magnitude. Power system oscillations become much worse 
following a large disturbance. In some cases the low frequency growing oscillations cause loss of power supply to custom 
loads, loss of synchronism among generators, or they reduce transmission capability of long transmission lines; these are 
expressions of power system instability. Voltage oscillations in a power system indicate a high degree of its vulnerability, 
because a change in the condition of power system could easily lead to a progressive drop in voltage at all buses in the 
transmission network, and to voltage collapse. Distribution networks are practically exempted from this problem due to their 
passive nature. Over the last three decades, the problems of low frequency oscillations in power systems have assumed 
importance. 

The stability issue of power systems can thus be stated in three aspects [3], [4]: (i) rotor (or power) angle stability, (ii) bus 
voltage stability and (iii) frequency stability. An upset in the balance between power generation and power demand can affect 
frequency stability; the frequency could no more be maintained within the stability limits. The inability of the power system 
to meet the demand for reactive power causes voltage instability. However, voltage instability does not always occur alone; 
often angle and voltage instabilities are associated. One may lead to the other and the distinction may not be clear. 
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Moreover, as power systems are nonlinear, their stability depends on both the initial conditions and the magnitude of a 
disturbance. Consequently, rotor angle and voltage stability can be divided into steady-state or small-signal stability and 
transient or large disturbance stability [1], [5], A distinction between these two types of stability is important for 
understanding the underlying causes of the problems in order to develop appropriate controls and operating procedures. 
Under steady-state operating conditions, bus voltages and line currents of all the synchronous generators must have the same 
frequency and the rotor mechanical speed of each generator must be synchronized to this frequency: the generators are said 
to operate in synchronism. A perturbation on the power system that results in a change of the electrical torque could cause 
oscillations of the generator rotor around its equilibrium point due to insufficient damping torque component, or a steady 
increase in the rotor angle due to insufficient synchronizing torque component. A high magnitude of rotor oscillations or the 
angle drift of the generator rotor can entail loss of synchronism with the rest of the network and disconnection of certain 
generators; rotor angle instability is associated with slow loss of synchronism among generators. 

In the Cameroonian power system, considering the southern interconnected grid (SIG) consisting of two hydropower plants 
with 16 turbine-generator units and one gas plant, which supplies six regions of the country where almost 90% of economic 
and industrial activity in the country is concentrated, it constitutes the core of the Cameroonian power system, and provides 
the bulk of Cameroons’ electricity [6]. In the SIG, power generation is centralized and hydro -dominant; the three power 
generating stations are connected to a very high voltage transmission substation in Mangombe, and from there power is 
transmitted to the different load centers. In the past two decades the SIG suffered from instability which prevented the grid 
from beeing fully utilized, and from very poor power quality; it experienced many perturbations due to power imbalance - 
insufficient power generation and increased active or reactive power demand. Low frequency oscillations were observed in 
bus voltages and line currents; they often caused loss of power supply to custom loads in the nearby region, or reduced power 
transmission capability of long transmission lines like those from Mangombe to Yaounde and Bafoussam. The transmission 
network in the SIG also often experienced voltage drops and voltage oscillations during periods of high power demand. 
During dry seasons, perturbations were worsened due to insufficient water flow for production; power furniture could not 
meet demand and was frequently interrupted for hours or days - load shedding was often used to match power demand with 
power production. 

In power systems, a particular issue encountered at the generating plant level is to maintain stability under various operating 
conditions [3]. The problem of small-signal stability in generating plants is usually one of insufficient damping of their 
oscillations. The first solution to this problem was to realize damper windings on generators and turbines [7]. But with 
generators equipped with slow excitation systems (DC and AC excitations), the available synchronizing torque was weak. 
This weakness of the synchronizing torque was found to be causing system instability when power systems started operating 
close to their stability limits. To enhance the synchronizing torque, fast excitation systems (static excitations) were used [5], 
To avoid this second problem, an automatic voltage regulator (AVR) was installed on the generators, so that it acted on their 
excitation system. The generator control is one of the basic control means of power systems; it incorporates also a speed 
governor for speed control. The AVR introduces a positive synchronizing torque component in excitation systems, and a 
negative damping torque component that reduces the damping of the system oscillations [3]. A poorly damped power system 
could be subjected to cascade failure after a variation of demand. The low frequency oscillations in power plants 
incorporating generators equipped with a static excitation system were traced to fast voltage regulation. The generators of the 
Song-Loulou power plant are equipped with static excitation systems; thus there is a risk that by greater amounts of power 
flowing across High Voltage or Very High Voltage long transmission lines an increasingly nefarious oscillatory instability 
factor be added compared to earlier situations, leading to low frequency oscillations in the SIG. 

To enhance the stability of power systems. Power System Stabilizers (PSSs) and Flexible Alternating Current Transmission 
Systems (FACTS) are being currently used as additional controllers. Although FACTS are quite effective on damping low 
frequency oscillations, they are merely used to enhance the capability of transmission lines, but they are not used in the SIG. 
The PSS is added to the AVR in the generator excitation system to enhance the damping of low frequency small-signal 
oscillations. The operation of the Cameroonian SIG equipped with classical AVR and PSS has been shown to be 
unsatisfactorily. Experience has shown that a classical power system stabilizer tuned for an operating point do not show good 
performance by perturbations in different operating points [8]. There is need for additional generator controllers which 
simultaneously increase the damping of power system oscillations and the synchronizing torque, in order to enhance the 
overall steady-state and dynamic stability (angle, voltage, and frequency) of the Cameroonian SIG. 
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Modem control theory with multivariate approach could offer some tools for this task. Until now, from the existing generator 
control systems none is the implementation of modern control theory [9], they were designed using frequency or root locus 
method, and they were designed to solve a single problem. It was seen that they helped solve one aspect of the problem, but 
to the detriment of another aspect. Control system design generally begins with the realization of an appropriate model for 
the dynamic system, which will be used to evaluate its overall performances. Power system operators and planners need 
effective tools to assess the overall stability of the system and take appropriate actions to enhance its stability over a wide 
range of operating conditions: the model simulation. Modeling of power plants and their loads is therefore fundamental to the 
evaluation of their performance by simulation and design of their control. 

Numerous research works have been done and published on modeling and simulation of power systems for small-signal 
stability studies. Demello and Concordia [4], Kundur et al. [11], and Chow et al. [12] developed a linearized model of 
synchronous generator and its excitation system connected to an infinite bus in the form of block diagram, in frequency 
domain. Baghani and Koochaki [8] developed a linearized state-space model of a synchronous generator and its excitation 
system. Their model was derived from the steady-state harmonic operation which is linear. Pedro Camilo de Oliveira e Silva 
et al. [13] developed a linearized state-space model of a complete power plant. 

For the Cameroonian SIG, state-space models of water flow system, turbine -generator unit, voltage control system, speed 
control system, and load connected at the generator terminals will be developed. Thank to the architecture of the SIG that is 
rather a radial than a meshed one, the SIG can be considered as simple, and will be modeled as a one area single machine 
system. The model must overcome some of the shortcomings of the aforementioned models. Thus, the turbine -generator unit 
model will be built by physical principles and put in state-space form. Small disturbances hypothesis permits the linearization 
of the system of equations and the application of all the wealth of the modern control theory. The other devices will be 
transformed from frequency representation form that is already linear, into state -space representation form. The performance 
of the integrated overall system will be assessed using MATLAB simulation. 

The rest of this paper is structured as follows: in section II the feedback control structure of the system is presented. In 
section III the detailed development of a linearized state-space model of the Cameroon SIG is realized; this linearized form is 
suitable for practical study of the classical control systems performance. Some illustrative simulation results are presented 
and discussed in section IV. In the final section of the paper some concluding remarks are mentioned. 

II. Feedback Control Structure of the Power System 

The Cameroonian SIG incorporates in reality many autonomous feedback control systems, each consisting of a generating 
unit (turbine and alternator), and the load. It can be modeled as a one area single machine system. 



Fig. 1: Feedback Control Structure of the power plant 


The structure of the feedback control system of the SIG is shown in fig. 1 [14]. A generating unit incorporates three feedback 
control loops for the regulation of the three main variables influencing the power output. These variables include the flow 
rate of water impacting the turbine, the alternator rotor speed and the alternator armature voltage which determines the 
magnetic flux permeating the armature circuit. Each of these variables is maintained at a reference value by an appropriate 
feedback control loop. 
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III. State-space Modeling of the Power System 


3.1 The Turbine- Alternator Sub-system 


From a system dynamics perspective, the alternator is the most complex component of the power system, due to the presence 
of several coils in armature and in rotor, which are mounted on separate geometrical axis and are magnetically mutually 
coupled. The alternators used in hydro electric power generating plants are salient poles synchronous machines. Their 
permeance varies across the air gap with respect to the pole axis (d-axis) of the polar wheel and to an axis at right angle to the 
pole axis, the quadrature axis (q-axis) [14], As a logical consequence of this, an orthogonal system with d-q axis is 
introduced on the rotor, with the d-axis superimposed on the polar axis and the q-axis lagging the d-axis. The various 
inductances in the alternator are represented by six coils [3], [15]. The stator contains three coils mounted on three axis 
separated from each other by an angle of 120°. The inductor coil F is located on the rotor in the d-axis. The damping effect is 
represented by the short-circuited fictitious coils KD and KQ in the d- and q-axis respectively. Some inductances vary with 
the angular position of rotor with respect to the axis of armature reference coil 9 which in turn varies with time. This 
introduces considerable complexity in the machine differential equations whose coefficients are time variable. The 
elimination of the time dependency of the inductance terms is obtained by Park’s transformation [15]. Due to this 
transformation, the magnetizing effect of the three-phase armature currents is replaced by that of their components i d (t) and 
i q (t) flowing in two imaginary coils located on the d- and q-axis of the rotor, and the magnetizing effect of the damping coils 
is due to the currents i KD (t) and i K Q(t) flowing in coils KD and KQ. 

3.1.1 Electrical State-Space Model of the Alternator 


The electrical state-space model of the three phase alternator in d-q-0 components is [14]: 

f did(t) 'j 
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iq(t) 

i F (t) | + |Br 
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( 1 ) 


where the zero-sequence voltage u 0 and current i 0 are zero, since the alternator operation is assumed to be in equilibrium. The 
system matrices [a g I and [b g J are independent of time and are defined in Appendix 1. 


3.1.2 Mechanical model of the Rotor of Alternator-Turbine Unit 


Assuming rigid coupling of turbine and alternator, the motion of the rotor is described by the following equation: 


Jt— f = T m -T em (2) 

dt 2 

where, 

J, = moment of inertia of all rotating elements 
T m = mechanical torque produced by the turbine 
T em = electrical torque opposing the motion of the rotor 
0 m = mechanical angle 

In the Park transformed electrical model of the generator, the electromagnetic torque results from the interaction between the 
total flux linkages and 'P q with the fictitious windings in d- and q-axis, and the currents through these windings i d and i q . 
It is expressed as follows [14], [15]: 
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Ten^P-Kl-iq-^q-id) 


(3) 


where, 

Yd and 'F (| = fluxes permeating the coils d and q, respectively 
i j and iq = currents in the coils 
p = number of pairs of poles 

The expressions of 'F d and 'F q are given by [14], [15]: 

= L d i d (t) + -J—Mps ipCt) + ^~ m KD , S i KD (t) 

< 

'PqOO = LqiqCOT J^M K Q S i K Q(t) 

The electromagnetic torque then has the detailed form: 

T em = pj( L d _L q )id iq + J~ M FS iq *F j + pj-^“ M KD,S iq if® ~J^ M KQ,S id ^ICQ 


(4) 


(5) 


Even at synchronous speed, the rotor oscillates and EMFs are then induced in the damping coils and in the rotor iron. The 
currents driven by the induced EMFs (eddy currents in the rotor iron) interact with the flux linkages of the damping coils and 
produce a torque which dampens the oscillations (Lenz’s law). The effects of eddy-currents in the rotor iron and of currents 
in the damper coils can be lumped into a damping torque. The second term in equation (5) is obviously the damping torque 
T d which is moreover proportional to the rate of change of the angle of polar wheel. It is described by the equation: 


T D -P]i “ M KD,S iq iffD +i “ M KQ,S M ^KQ f- k D 


d8p 

dt 


(6) 


Substituting for T em and T D in equation (2) results in the non-linear differential equation that describes the mechanical 
behavior of the generator rotor and the turbine: 


d8 p | pk D d5 p 
dt J t dt 


P- T _P_( L 
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d 2 0 d 2 8 p d5 p 

where = = . 

dt 2 dt 2 dt 

3.1.3 Linearization of the mechanical model 


During stable and balanced operation of the alternator, the stator phase voltages are sinusoidal of constant magnitude V , 
they form a balanced three phase system. Applying the inverse Park transform to the system of phase voltages results in the 
components [14] : 


fu d (t) = — Vj V sin 8 p = -U sin 8 p 
juq(t) = a/ 3V COs8p = U COsSp 


( 8 ) 


When the system oscillates about an operating point (u d0 , u q0 , Uo, 5 p o) with the small variations AU and A5 p of the variables U 
and 5 P respectively, the small variations of the components u d and u q can be approximated as: 


Au d « -sin S p o AU - (Uq cos 8 p o ) A8 f 
Au„ ®cosS„o AU - (lJo sin 8 p0 )A8p 


(9) 
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Linearising equation fl) and ignoring the damping currents i KD and i K Q whose effect has already been taken into account by 
introduction of the damping torque T D , results in the following state equations of the three phase alternator: 

f dAi r ' 


dt 

dAi q 

dt 

dAip 

dt 
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b g ]- 
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(10) 


where the system matrices [Aq] and [Bq ] are given in appendix 1. Equation (7) can also be linearized by considering 
small variations about an operating point (i d0 , i q o. h o, n 0 , T m0 , P m o)- The linearized version of this equation is of the form : 
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(11) 


J 


Combining the state equations (10) and (11), the state-space representation of the turbine -alternator sub-system is obtained in 
the form: 
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( 12 ) 


where the system matrices [A jq and Bjq are given in appendix 1 . 


3.1.4 Model of the Hydraulic Turbine 

The water gate opening Y and the head H t have a direct influence on the water velocity in the penstock and thus on the flow 
rate. In stability studies of power systems, if mechanical power P m is the output variable and gate opening Y the input 
variable, the non ideal hydraulic turbine can be represented by a transfer function F TH as follows [3], [15] : 


F XH (s) = K t 


l + b thl T w s 

l + a thl T w s 


(13) 


where, 

s= Laplace operator 

T w = hydraulic time constant for a given load level 

K t = gain of hydraulic turbine 

a th i and b thl = coefficients of transfer function 

The canonical state-space representation for small variations about a stable operating point, derived from the function (13) is 
of the form: 
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where the coefficients of the state-space model are represented in Appendix 2. 

3.2 State-Space Models of the Control Systems 

As shown in fig. 1, there are three control loops in the system -the flow, the speed and the voltage control systems. 

3.2.1 The Flow Control System 

The block diagram of the flow control system [16] is shown in fig. 2. 



Fig. 2 : Block Diagram of the Flow Control System 

The system incorporates a PID controller. The variables and parameters in the block diagram are defined as follows: 

Q = water flow rate 

Q r ef = reference water flow rate 

K P = proportional gain of the PID flow controller 

Tj = rate time of the PID flow controller 

T, = reset time of the PID flow controller 

T f = filter time constant of the flow-head converter 

K b = rate of slope 

The closed-loop transfer function of the system is of the form: 
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(15) 


where the coefficients a j rh and h irh are represented in Appendix 3. The corresponding canonical controllable state-space 
representation is of the form: 
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(16) 
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3.2.2 The Speed Control System 

The block diagram of the speed control system is shown in fig. 3 [3], [17]. The system incorporates a PI controller. 



Fig. 3: Speed Control System 

The variables and parameters in the block diagram are defined as follows: 

Y : W ater gate opening 

n: Angular speed of the alternator 

P: Active power supplied to the load by the alternator 

n ra : Speed set point 

P r( sf : Active power set point 

T n : Time constant 

K c : Proportional gain of PI 

Kj : Integral gain of PI Controller 

B[ : gain of limiter 

T x : Time constant 

K s : Gain of the servo valve 

T s and T v : Time constants of the servo valve and valve-wagon units respectively 

B t : Transient droop 

R p : Permanent droop 

Tii : Minor loop Time constant 

In addition to the main speed control feedback loop, a subsidiary power feedback loop is applied to the system to limit the 
active power supplied during transient operation or in situations of overload. 

The closed-loop transfer function of the system is of the form: 


b 0rv + b lrvn s + b 2rvn s 

a 0rv + a lrv s + a 2rv s2 + a 3rv s3 + s4 
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+ b 
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(Pref ~P) 


(17) 


where the coefficients a lrv , bi rvn and b irvp are represented in Appendix 4. The canonical state-space representation, derived 
from the closed-loop transfer function, is of the form: 
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3.2.3 The Voltage Regulator 

The block diagram of the voltage regulator is shown in fig. 4. 


u 


Fig. 4 : Block Diagram of the Voltage Control loop 

The parameters in the block diagram are defined as follows: 

Uref : Armature voltage set point 
U : RMS value of armature voltages 
I F : Excitation current 
U F : Excitation voltage 

U r : Voltage representing the RMS value of Armature coils voltages 

Ki : Gain of the voltage regulator 

ii : Time constant of the voltage regulator 

K 2 : Gain of the reference voltage amplifier 

K 3 : Slope of the thyristor characteristic 

K 5 : Gain of the excitation voltage regulator 

if : Time constant of the excitation voltage 

K 6 : Gain of the current regulator 

i6 : Time Constant of the current regulator 
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K 8 : Conductance of the armature voltage control chain 
is : Time constant of the armature voltage control chain 


The canonical state-space representation. 
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where the coefficients a lrt and b rt are represented in Appendix 5. 

3.2.4 The Power System Load 


(19) 


The behavior of the power system load has great influence on its stability and that of the generating units. The multiplicity of 
loads and the diversity of their behavior make their modeling a difficult task. Traditionally, load models are classified in two 
broad categories: static and dynamic loads. Here we consider static loads whose voltage-frequency behavior can be described 
by exponential equations [3], [15] of the forms: 


P L = P L o 


' u v " 


i u o; 


(l + KpfAf) 


Ql = Qlo 


a u y [ \ 

TT (l + KqfAf) 

V u oJ 


(20) 


where : 

P L0 = initial value of active power in pu 

Q lo = initial value of reactive power in pu 

U 0 = initial RMS value of armature voltages in pu 

U = RMS value of generator armature voltages in pu 

Af = variation of the frequency in pu 

m, n = voltage exponent for active and reactive loads respectively 
K pf = frequency sensitivity for active load 
K qf = frequency sensitivity for reactive load 

The linearized version of equations (20) is then: 

mP 

AP L =- r 7^AU+P L0 K pf Af 

u 0 

AQ L =^ MI AU + Q L0 K qf Af 
u 0 


( 21 ) 


The small variation of the RMS value of the armature voltages in the reference frame of the alternator is of the form : 
AU = R d Ai d + R q Ai q + U 5 A8 p + U f Af 


( 22 ) 


where the coefficients R d , R q , Ug and Uf are represented in Appendix 6. 

3.3 State-Space Model of the Power System 

Combining the models of the constituent sub-systems, we obtain an 18 th order state-space model of the power system that is 
of the form: 
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(23) 


where, 

(Ax) = 18 x 1 state vector 

(Ay) = 13x1 output vector 

[Ag] = 18x18 matrix (the A-matrix) 

[B g ] = 18x3 matrix (the B-matrix) 

[C S ] = 13x18 matrix (the C-matrix) 

The elements of these different matrices and vectors are presented in Appendix 7. 

IV. Simulation Results and Discussion 

In this paper, to show the validity of the linearized state-space model of the Cameroonian SIG, and hence to investigate the 
performance of its PID controller, PID controlled voltage and speed regulators, simulation works have been carried out for 
the single machine system connected to a load, using MATLAB programs designed by the authors. Most of the model 
parameters used in the ensuing simulation are given in the appendix and were obtained from technical manuals of the system. 
Some parameters were obtained from experiments on a prototypal alternator [18], 

In order to prove the validity of the model, the results are compared with those given in source books like [3] and [15], The 
perturbations considered in this paper are input unit steps which are applied to the system. We consider the operating point 
given in appendix 8 which also gives the values of sub-systems parameters. The parameterized state-space model is 
presented in Appendix 7. The simulation results are step responses of water flow deviation, armature voltage deviation, rotor 
speed deviation and rotor angle deviation. They are depicted respectively in figures 5 to 8. 

Fig. 5 illustrates the step response of the dam water flow system given in fig. 2, to a unit step disturbance of the reference 
water flow rate. 



Temps (sec) 

Fig. 5 : Response to the set point unit step of the water flow 

The PID controller really stabilizes the system. The dynamics of the water flow rate deviation exhibits during post 
perturbation state many oscillations with large initial amplitude which then decreases slowly, before it settles to its steady 
state value 1,0. The system is poorly damped and it is not regulated back to zero deviation after a unit step perturbation of the 
water flow set point. The system with classical PID controller has a long settling time and a bad accuracy. 

Fig. 6 illustrates the step responses of the turbine -alternator system, following a disturbance of the reference armature 
voltage. It shows the deviations of the armature voltage AU, the rotor speed An and the polar wheel angle A5 P . 
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Temps (sec) Temps (sec) Temps (sec) 


Fig. 6: Responses to the set point unit step of the armature voltage 

The dynamics of the armature voltage deviation exhibits a first swing with overshoot during post perturbation state, before it 
decreases, oscillates slowly and settles to its non zero steady state value; it contains high frequencies due to transients in 
stator voltages. The PID controlled AVR really stabilizes the system, but not accurately; the armature voltage perturbation is 
poorly damped and not regulated back to zero by the AVR. The dynamics of the rotor speed deviation displays during post 
perturbation state, an acceleration of the rotor in the first half cycle and its deceleration in the second half cycle and so on, 
with decreasing amplitude; it oscillates slowly and settles to its zero steady state value; the perturbation is completely cleared 
by the speed controller. The dynamics of the rotor angle deviation exhibits during post perturbation state a first positive 
swing followed by a negative swing, before it oscillates slowly after a negative overshoot and settles to its non zero steady 
state value. 

Fig. 7 illustrates the step responses of the turbine-alternator system, following a disturbance of the reference rotor speed. It 
shows the deviations of the rotor speed An, the rotor angle A5 p and the armature voltage AU. 



Temps (sec) Temps (sec) 


Fig. 7 : Responses to the set point unit step of the rotor speed 

The dynamics of the rotor speed deviation exhibits during post perturbation state, first a decrease followed by an increase and 
oscillations with decreasing amplitude, before it settles to its zero steady state value; the perturbation is completely cleared 
by the speed controller. The rotor is decelerated in the first half cycle and accelerated in the second half cycle and so on, with 
decreasing amplitude; it oscillates slowly and settles to its zero steady state value. The PID controlled speed regulator really 
stabilizes the system, accurately. The rotor angle deviation first exhibits a decrease and then it increases steadily to reach a 
positive overshoot; it then decreases and oscillates to stabilize at a non zero value. The armature voltage deviation first 
exhibits a decrease and then it increases steadily to reach a positive peak value and decreases to stabilize at a non zero value. 
The PID controlled AVR really stabilizes the system, but not accurately. The initial retardation of the rotor has an aftereffect 
on the rotor angle and the armature voltage. They are caused by high frequency transients in the transformer voltage terms in 
the stator voltages. 

Fig. 8 illustrates the step responses of the turbine -alternator system, following a disturbance of the reference active power 
demand. It shows the deviations of the rotor speed An, the polar wheel angle A5 p and the RMS value of the armature voltages 
AU. 
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Fig. 8 : Responses to the set point unit step of the active power demand 

The dynamics of the rotor speed deviation exhibits during post perturbation state a high negative peak that means a sharp 
deceleration in the first half cycle, followed by an acceleration in the second half cycle and so on; it increases to a high 
positive overshoot and then decreases and oscillates around zero, before it settles to its zero steady state value; the power 
perturbation is completely cleared by the speed controller: the PID controlled speed regulator really stabilizes the system, 
accurately. The dynamics of the rotor angle deviation and the armature voltage deviation exhibit a first decrease to negative 
values (backswing phenomenon) and then a steady increase to positive values. They exhibit overshoots and stabilize at non 
zero values. The PID controlled regulator does not clear the perturbation of the power demand. 

The high frequency oscillations affecting the rotor speed are due to the derivative action of the rotating mechanical system of 
the turbine-alternator system. 


V. Conclusion 

A linearized dynamic model of the Cameroonian southern grid has been developed, using the state-space approach. The 
model is modular in structure. Intensive MATLAB simulation has been carried out to investigate the behavior of the southern 
grid of the Cameroon power system, after small disturbances. The simulation results have shown that classical controllers 
actually in the southern grid of the Cameroon power system ensure stability of the SIG, but they are not accurate, they are 
slow. The power system is poorly damped. This assessment can serve as a basis for the design of most efficient additional 
automatic controls to operate in combination with classical controllers present in the SIG. 
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APPENDIX 

Appendix 1: Matrices of Generator and Turbine-Generator Unit Models. 

[ A G J = [ L G ] 1 '[ R g] an d [b g J= [ L g] 

where, 



L d = inductance of the imaginary coil in the d-axis 
L q = inductance of the imaginary coil in the q-axis 
L F = inductance of the inductor coil F 
Lkd = inductance of the damping coil KD 
L kq = inductance of the damping coil KQ 

M FS = maximum mutual inductance between the inductor coil F and the stator coils 
Mkd.s = maximum mutual inductance between the damping coil KD and the stator coils 
M K q,s ^maximum mutual inductance between the damping coil KQ and the stator coils 
M Fi kd = maximum mutual inductance between the inductor coil F and the damping coil KD 



L F sin5p 0 

L f U 0 cos5 p0 

Ldh-^FS 

L d L F -^ M ?s 

- cos8p 0 

U 0 sin5 p0 

L q 

L q 

- J-Mfs cosS p0 

_ h MfsU ° s ' n ^ p ° 

L d L F -^l\4 

L d L F “|M^ s 
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Appendix 2: Parameters of Hydraulic Turbine. 

Non ideal Hydraulic Turbine: K t =1,5 ; a thl3 = 1,1 ; 

a t = 2,064 • 1CT 10 x P L0 + 0,57 ; a th21 =4,541 -lCT 9 xP L0 +1,18 ; b t =a t _ a thi3 x a th 2 i 
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Appendix 3: Coefficients of Transfer Function of Water Flow Control System. 
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Appendix 4: Coefficients of Transfer Function of Rotor Speed Control System. 
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Appendix 5: Coefficients of Transfer Function of Voltage Control System. 
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Appendix 6: Coefficients of Power System Load Model. 
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The voltages u D0 and u Q0 , and the currents i D0 and iQ 0 , are Kron’s components of armature voltages v a , v b and v c , 
and currents i a , i b and i c respectively, in a new coordinate system D, Q and 0. They represent steady-state voltages 
across and currents in two imaginary coils located on the D and Q axis of the stator, for a given operating point, 


and are obtained as follows: 
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The Kron’s transformation matrix [K(0 O )] is similar to the Park’s. 
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Appendix 7: Parameterized State-space Model and Elements of Model Matrices. 


A 1,1 -LpCQb 
A 1,4 =L F co b 

A 2,1 =co b 


— R s + R d sin Yo 
L d L F -M 2 s 
Up cos Yo +U 5 s in Yo 
L d L F -M 2 S 
co sO L d- R q cos Yo 


A 1,2 “ L F®b 
A 15 = L F co b 


-ro s 0 L p + R q siny 0 


L d L F -M^s 
pU F sin y Q 
L d L F -M^s 


Ai 3 - R F co b 


\/3/2 -M 


FS 


A 2,2 - ro b 


-R s — Rq cosy 0 


A 115 __ V3/2Mp S co b 


L d L F -Mfs 

b 0RT 
L d L F ~ M FS 


A 2 ,3 - © b 


a /3/2M Fs (o s o 


A 2 4 - K> b 


U 0 siny 0 -U 5 cosy 0 


a „ pUp cos Yo A . /^Yaa R s - R d cosy 0 

A 2,5- _ ®b ; A3. -co b V 3 /2Mp S — 


LdLp-Mfe 


A„- Mb -^7lM -at- 


RpL 


L d L F -M^s 


F L d 


L d L F -M^ s 


Ay.-avAAM Fs PHi£5!l!> 


LdLf-Mps 


L d L F -M^s 


^3,15 “ co b 


L ri b 


d u 0 rt 


L d L F -Mp S 


A 4 5 - co b 


A 5,l - _K i 


A 5,2 “ _K iq 


A 5,3 “ _X iF 


'7,8 


'10,4 


= -K n 

A 5,6 = 

: K pm b 0TH 

L =K pm d TH 

b 02 rv 

A 5,12 

= koirv 

A 6,8 

= b llrv 

= 1 

A 8,9 

= 1 


A 5,7 - K pm d TH b 01rv 


A 5,8 - K pm d TH b llrv 


A 5,13 “ K pm d TH b 22rv 


A 6 ,ll -bo 2 rv 


A 6,12 -bl 2 rv 


'9,10 


= 1 


mP m R 


' 10,1 


LO r 'd 


noPioUg 


'10,7 - a 0rv 


' 10,8 “ a lrv 


U 0 

A 


A 6,6 — a OTH 


A 6,13 — ^22rv 


' 10,2 


mP L0 R q 


10,9 “ a 2rv 


' 10,10 “ a 3rv 




International Journal of Engineering Research & Science (1JOER) 


ISSN, [2395-6992] 


[Vol-2, Issue-8, August- 2016] 


11,12 “ 1 

A 12,13 = 1 

A 13,14 = 1 

A 14,5 = 

-1 A 14,11 = — a 0rv 

A 14,12 = _a lrv 

14,13 = _a 2rv 

A 14,14 = _a 3rv 

A 15,15 = 

a Ort 

A 15,16 = _a lrt 

A 15,17 = _a 2rt 

16.3 = 1 

18.4 = LI 5 

A 16,16 = a 3rt 
A 18,5 =PUf 

A 17,17 = _a 4rt 
A 18,18 = a 6rt 

A 17,18 

= — a 5rt A 18,l =R d 

A 18,2 = R q 


co s0 = 2 nf, 


N > 


Yo = atari 


' 0 soLql 0 cos(p-R s l 0 sirup 

U 0 +R S I 0 cos(p + co S 0 L q l 0 sirup 


B io,i - 1 


B 14,2 - 1 


B 17,3 


= 1 


Ci,i -1 c 2,2-! C 3;1 -R Dd cosy 0 +X Qd siny 0 c 3,2 - x Dq cos Yo _R Qq sin Yo 

c 3,4 =U Dd cosy 0 +U Qc |Siny 0 -u q0 C 3<5 = p(-U Df cosy 0 -U Qf siny 0 ) 


c 4,l = R Dd sin Yo - X Qd cos Yo C 4 2 =X Dq siny 0 +R Qq cosy 0 

c 4,4 =U Dd siny 0 -U Qd cosy 0 -u q0 C 4<5 =p(-U Df siny 0 +U Qf cosy 0 ) 


C 5,3 - 1 


C 6,15 -^Ort 


C 7 ,l “ R d 


C 7,2 “ R q 


C 7 , 4 =U S 


C 7,5 -P U f 


_ 180 _ 60 

L 8,4 - '- 9,5 — ®b — r— 

n p- 2 n 


C 10,5 -P'fb 


C ll,l-“( L d L q) 


q ' qo 


' 11,2 


' 12,4 


( L d Lq) 


3 

. mP L o u 6 


U 


0 



- + P L 0 K pf 


^113 — Mpci 


mP m R 


1 FS 1 qO 


' 12,1 


LO r 'd 


Ur 


' 12,2 


mP L0 R q 


n Q L o R d 


' 13,1 


n Q L0 R c 


U r 


- 13,2 


' 13,4 


nQ L o u 5 

Un 


( 


C 13,5 -P 


n Q|_oUf 

v u 0 


+ Q10K 


LO l 'qf 


C 14,7 -L) 0 irv C 14 , 8 -' 3 llrv C 14 ,ll “ L> 0 2 rv C 14,12 “ L) 12rv C 14,13“L 


a l,l 

a l,2 

a l,3 

a l,4 

a l,5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

a l,15 

0 

0 

0 

a 2,l 

a 2,2 

a 2,3 

a 2,4 

a 2,5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

a 3,l 

a 3,2 

a 3,3 

a 3,4 

a 3,5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

a 3,15 

0 

0 

0 

0 

0 

0 

0 

a 4,5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

a 5,l 

a 5,2 

a 5,3 

0 

a 5,5 

a 5,6 

a 5,7 

a 5,8 

0 

0 

a 5,ll 

a 5,12 

a 5,13 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

a 6,6 

a 6,7 

a 6,8 

0 

0 

a 6,ll 

a 6,12 

a 6,13 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

a 10,l 

a 10,2 

0 

a 10,4 

a 10,5 

0 

a 10,7 

a 10,8 

a 10,9 

a 10,10 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

-1 

0 

0 

0 

0 

0 

a 14,ll 

a 14,12 

a 14,13 

a 14,14 

0 

0 

0 

0 

a 15,l 

a 15,2 

0 

a 15,4 

a 15,5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

a 15,15 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

a 16,16 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

a 17,17 

a 17,18 

a 18,l 

a 18,2 

0 

a 18,4 

a 18,5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

a 18,18 


22 rv 




International Journal of Engineering Research & Science (1JOER) 


ISSN, [2395-6992] 


[Vol-2, Issue-8, August- 2016] 


1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

c 3,l 

c 3,2 

0 

c 3,4 

C 3,S 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

c 4,l 

c 4,2 

0 

C 4,4 

c 4,5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

b rt 

0 

0 

0 

C 7,l 

c 7,2 

0 

c 7,4 

c 7,5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

c 9,5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

c 10,l 

c 10,2 

c 10,3 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

c ll,l 

c ll,2 

0 

C ll,4 

c ll,5 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

c 12,l 

c 12,2 

0 

c 12,4 

C 12,S 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

c 13,7 

c 13,8 

0 

0 

c 13,ll 

c 13,12 

c 13,13 

0 

0 

0 

0 

0 


r Ai d ^ 


'0 

0 

0~ 


Ai q 


0 

0 

0 


Aip 


0 

0 

0 


( Ai d ^ 

ASp 


0 

0 

0 



An 


0 

0 

0 


Aiq 

x t 


0 

0 

0 


Au d 

x rvl 


0 

0 

0 


AUq 

x rv2 


0 

0 

0 


Ai F 

x rv3 

[B S ] = 

0 

0 

0 


Au f 

x rv4 


1 

0 

0 

(Ay) = 

AU 

x rvl 


0 

0 

0 


A5 n 

x rv2 


0 

0 

0 


P 

x rv3 


0 

0 

0 


An 

x rv4 


0 

1 

0 


AT em 

x rtl 


0 

0 

0 


AP l 

x rt2 


0 

0 

0 


aq l 

x rt3 


0 

0 

1 


v A y 2 

V x rt4 ) 


0 

0 

0 



(Au) = 

AP r ef 1 

An r gf 


v AU r efJ 


Appendix 8: Values of Parameters of subsystems of State-Space Model. 

System Operating Point: 

U 0 = 1.2247p.u. ; P LO =0.9p.u. ; Q L0 =0.5578p.u. ; cos(p 0 =0.85 ; l 0 = 0.7487p.u. ; 

i F0 =1.094p.u. ; ij G0 = 0.97878 ; n 0 = 120tr/min. ; 8 p0 =21° ; T w =4.9661sec. 

Synchronous Generator: 

L d = 1.4938p.u. ; L q = 0.9959p.u. ; L CT =0.167p.u. ; M FS = 1.3268p.u. ; p = 25 

L f =1.852p.u. ; R s = 0.0152p.u. ; R F = 0.0002654p.u. , J t = 8800000kgm 2 . 

Water Flow Control System parameters: 

K b = 50 ; T f = 0,916 ; K p =l ; T d =5 ; T,=10 ; 

Speed Regulator System parameters: 

n ref =50 ; T n =2,22 ; K c =2,08 ; K; =0,416 ; R p =1,0 ; B, =0,44 ; T x = 1,344 ; K s =0,2 ; T d =2,02 ; B t =0,55 ; 
T s =0,07 ; T v =0,2 . 

Voltage Regulator System parameters: 
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K-l =0,1 ; x 1 =0,04 ; K 2 =5,39 ; K 3 =10,47 ; K s =1,486 ; i f =0,93 ; K 6 =0,971; x 6 =0,01 ; 
K 8 =0,000971 ; x 8 =0,01 . 

Power System Load: 

m=2 ; n = 2 ; K pf =3 ; K qf =-2. 




